Balls-in-boxes condensation on networks 
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We discuss two different regimes of condensate formation in zero-range processes on networks: 
on a g-regular network, where the condensate is formed as a result of a spontaneous symmetry 
breaking, and on an irregular network, where the symmetry of the partition function is explicitly 
broken. In the latter case we consider a minimal irregularity of the q-regular network introduced 
by a single Q-node with degree Q > q. The statics and dynamics of the condensation depends on 
the parameter a = \nQ/q, which controls the exponential fall-off of the distribution of particles 
on regular nodes and the typical time scale for melting of the condensate on the Q-node which 
increases exponentially with the system size N. This behavior is different than that on a (/-regular 
network where a = and where the condensation results from the spontaneous symmetry breaking 
of the partition function, which is invariant under a permutation of particle occupation numbers on 
the q-nodes of the network. In this case the typical time scale for condensate melting is known to 
increase typically as a power of the system size. 
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The formulation of the principles of non-equilibrium statistical mechanics is a challenge for theoretical physics. 
Non-equilibrium effects play an important role in many phenomena but we do not have a consistent theory which 
would describe them. An exception may be systems close to equilibrium for which we can gain some insight 
into their dynamics using ideas of the linear response and of the fluctuation-dissipation theorem. The far-from- 
cquilibrium dynamics is an uncharted territory. It is therefore useful to look for solvable models belonging to 
the latter class, which could teach us about what happens in this case. 

In this note we shall discuss zero-range processes on networks which belong to this class. On the one hand, 
they reveal an interesting, non-trivial and very rich behavior including a condensation, far-from-equilibrium 
dynamics and non-linear effects as for example the formation or melting of the condensate. On the other hand, 
due to a relation to the balls-in-boxes model [l|, [2| , which is analytically solvable, also these non-trivial effects 
are analytically treatable. 

The zero-range process is a stochastic process which describes a gas of identical particles hoping between 
neighboring sites of a lattice or network on which the particles reside (for reviews, see 0, 3 IE])- The transition 
rate for particles to hop from one site to an adjacent site depends only on the state of the node from which the 
particle hops and is independent either of the destination node or any other node. Since the hop rate requires 
only ultra-local information, the corresponding process is called zero-range process. 

When the density of particles exceeds a certain critical value, the system undergoes a condensation [l[ in 
which a condensate is formed on a single node of the network which attracts a large number of particles. The 
condensation takes place in real space and not in momentum space as for the Bose-Einstein condensation. This 
type of condensation is called balls-in-boxes condensation, in short B-in-B or backgammon condensation because 
it was first discovered in the balls-in-boxes (backgammon) model [j], 0] • 

The B-in-B condensation is observed in many systems. For example in statistical models of random trees 
(also called branched polymers) one observes a phase transition between a phase of generic elongated trees to 
a phase where a typical tree looks rather like a bush with a singular node which has a finite fraction of all 
branches [3, @| . The statistics of the degree distribution in the tree model can be mapped onto the B-in-B 
model. The emergence of a singular node on a tree corresponds to the emergence of a condensate of many 
balls in one box in the B-in-B model. A similar geometrical phase transition is observed in models of quantum 
manifolds discussed in the context of quantum gravity 0, Isf- At the phase transition the quantum manifolds 
collapse to a very singular geometry whose volume is almost fully concentrated in the closest neighborhood of a 
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single point of the manifold. More precisely, the ratio v/V of the volume v of the neighborhood within a radius 
of order of an ultraviolet cut-off around this point to the total volume V of the whole manifold is finite in the 
limit V — > oo even if the radius (ultraviolet cut-off) is kept constant. Again this phenomenon can be viewed as 
a particular realization of the B-in-B condensation Q • 

The B-in-B condensation explains many other phenomena as for instance wealth condensation Q , emergence 
of singular nodes in complex networks [l(J EJ , emergence of the Hagedorn fireball in hadron physics [12, EH 
and some transitions observed in shaken granular gases [14 , flSj ] or the phase transition in zero-range processes 
21]. Actually, the mathematics of the condensation is also almost identical to that of the Berlin-Kac phase 
transition in the spherical model [l6l |. In other words the B-in-B mechanism is quite generic and common. 

The state of the zero-range process is characterized by the distribution of the numbers of particles at all N 
nodes of the network: {rrii} = {mi, . . . , mjv} where rrii denotes the occupation number of the i-th node. The 
total number of particles M — m\ + . . . + is constant during the process. Particles hop from non-empty 
sites with a rate u{m,i) depending only on the site occupation number m,. The outgoing current of particles 
from a site i is distributed equally among all qi links emerging from the node, so the effective hop rate per link 
is u{m,i)/qi, where qi is called node degree. The function u(m) is identical for all nodes, but the factor \ jq\ is 
not since it explicitly depends on the node degree. In this note we assume that the network topology is fixed 
and so is the degree sequence {q{\. The most fundamental question is whether the process has a steady state 
and if so, whether it is unique. The answer to this question is affirmative if the network is connected. In this 
case the process has a unique steady state which depends only on the node degrees {qi} and the numbers of 
nodes N and particles M. This steady state corresponds to the only equilibrium state which is sooner or later 
reached by the process. In equilibrium, the partition function can be calculated analytically 
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where 
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The weight p(m) is identical for every node, but the total node weights Pi(m) = p(m)<7™ have an additional 
contribution q™ explicitly depending on the node degree. Because of the presence of the (Kronecker) delta 
function under the sum in Eq. ([T]) , the partition function does not entirely factorize into a product of independent 
weights for individual nodes. The constraint on the total number of particles plays an important role as we 
shall see below, because the occupation numbers of individual nodes are not independent of each other. 

The statics of the zero-range process is equivalent to a B-in-B model with the partition function (JT|) describing 
a system of M identical balls distributed in N boxes, each having a weight function pi(m). The probability 
that in equilibrium the system is in a state {m,} reads: 
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P(mi,...,m N ) = -^Ylpim-Mr' = ^Y[Pt( m t), (3) 
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where as before mi + • • • + mjsr = M. It is interesting to notice that this probability is invariant with respect 
to the following change of the weights: 

p(m) -» Ce" m p(m). (4) 

Indeed under the change ^ the partition function (fT]) changes as Z C N e> lM Z. The multiplicative factor 
which appears in front of the partition function cancels out in the probability ^ and thus the statistical 
averages do not depend on the parameters C and p.. This invariance is an important property of the model. 
The parameter C has the meaning of a normalization and can be used for example to normalize the weights 
to a probability. One should notice that the hoping rate u(m) = p(m — l)/p(m) is not affected by C. The 
parameter fi or more specifically e M rescales the hoping rate u(m) — > e _A1 u(m) or equivalently stretches the time 
scale t — > te^. 

The probability distribution (|3|) encodes the full information about the static properties of the system. For 
example one can calculate statistical averages of any observable X: 
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or correlations (XY) — (X)(Y) etc. A particularly interesting observable is the number of incidents that the 
i-th node is occupied by m particles: 

71* (m) = (S mum ) = Pi(m) Zt ^ (6) 

where, using the B-in-B analogy, Pi(m) is the weight of the i-th box, Z is the partition function for the total 
system of M balls in N boxes and Zi(m) is the partition function for M — m balls in the N— 1 remaining 
boxes. In a similar manner one can count multi-node distributions and average them over all configurations. 
For example averaging the two- node incident function S mitm 8 m . jn over all configurations gives the probability 
iTij(m,n) that in equilibrium there are simultaneously to particles at node i and n particles at node j. 

A zero-range process is said to be in the condensed phase if a finite fraction of all particles tends to occupy one 
node or, if one phrases it in terms of the underlying B-in-B model, if a finite fraction of all balls is concentrated 
in one box. The effect depends on the density of particles per node (or balls per box) p = M/N. If one keeps p 
constant and takes the limit N — > oo, one sees that above a certain critical density p c one box contains on the 
average a large number of balls (p — Pc)N which grows with N, while any other box has only p c balls. 

The most thoroughly studied and probably most surprising example of the B-in-B condensation is the con- 
densation which takes place in a system of identical boxes, that is for which all the weight functions in Eq. 
([21) are identical: pi(m) = ■ ■ ■ = p^{rri). The corresponding zero-range process is realized on a g-regular graph 
which has identical degrees of all nodes. In this case the partition function ((TJ) is symmetric with respect to the 
permutation of the box-occupation numbers {mi}. The condensation appears there as a result of a spontaneous 
symmetry breaking which selects one out of identical boxes for the location of the condensate. The criterion for 
the appearance of the condensation is controlled by the asymptotic value of the hop rate u(m) for to — > oo, 

=u[m)^u 00 . (7) 
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For Moo = oo the corresponding critical density is infinite [22| and there is no condensation in the model. The 
system is always in a fluid phase. This can be intuitively understood because if the hop rate u(m) — > oo, 
it amounts to an effective repulsion between particles which in effect avoid to occupy the same site. On the 
contrary, for = there is no price whatsoever for particles to pay for a numerous occupation of the same 
site, so in effect the particles tend to condense. The critical density is zero in this case and therefore the system 
can be then only in the condensed phase. One can say that due to an effective attraction, particles tend to 
keep as close as possible to each other. The most interesting case is when Uqo is finite: 0<Uoo <oo. Actually 
it is sufficient to consider only the case Moo = 1 since within this model all other values of Uao are equivalent to 
one as follows from the invariance with respect to the transformation (U). In the remaining part of the paper 
we shall therefore stick to the choice itoo = 1. The quantitative behavior depends on the exact form of u(m) 
but the critical properties, such as critical exponents, depend only on how u(m) approaches unity when m 
goes to infinity. When it behaves as u(m) — 1 + /3/m + . . . for large to, then p(m) behaves asymptotically as 
p(m) ~ to"' 3 . In particular, one can choose p(m) to be 

= W l )m ml ^ (J3 - l)T(P) 
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When the density p — M/N exceeds the critical value p c , a condensate with NAp particles is formed P. 
where Ap = p — p c . The critical density is given by the formula 
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and can be concisely expressed in terms of the generating function K(p) — £ m _p(TO,)e -A " Tl as p c = —K'(0)/K(0). 
In particular, for p(m) as in Eq. (JHJ), the critical density is p c = l/(/3— 2). An interesting choice of the weights 
p(m) is when one demands that every box has at least one particle, and p(m) — to - ' 3 . In this case the generating 
function has the following integral representation: 
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which uncovers mathematical similarities between the B-in-B and the Bose-Einstein condensation. The statics 
of the model can be solved analytically [H, 0|, with critical properties depending on /3. 



The dynamics of the condensation has also been studied analytically 0, [l?], Ell ■ Two questions can be posed: 
what is the typical time scale for building the condensate from a homogeneous distribution of particles and 
what is its average life time. Here we shall focus on the latter. Once the condensate is formed, it moves 
across the system. It spends a long time at one particular site but sometimes melts and is rebuilt at another 
site. A typical time scale for melting the condensate has been derived using mean-field arguments [T^. In the 
mean-field approach one monitors only a single node of the network and derives effective equations balancing 
the in-flow and out-flow of particles for this node. One does not care about what happens in the remaining 
part of the system which in this approximation is treated as a homogeneous reservoir of particles, where 
fluctuations are much faster than the dynamics of the condensate. The monitored node is characterized by the 
distribution of the number of particles 7Tj(m) which, in an adiabatic approximation, is assumed to be that of 
the steady state. For a homogeneous system, 7ii(m) is identical for all nodes, so it is equal to the average over 
all nodes: 7Tj(m) = Wlim) = 7r(m), and it does not matter which node is monitored. The full information is 
encoded in ir(m). The occupation number of the monitored node may change in one step by one unit or stay 
constant. This sequence of changes can be viewed as a random walk in the effective one-dimensional potential 
V(m) = — hi7r(m). The waiting time for a condensate to melt can be thus viewed as the time needed for a 
particle to randomly walk from m*, which corresponds to the value of the condensate to* = NAp, to some 
too <C to*. This time is related to going through the maximum of the potential, whose position corresponds 
to the position of the dip of the function ln7r(m). This position is very close to to*/2 because the excess of 
particles is shared mainly by two nodes |l7]. In the condensed phase the shape of the function tt(to) can be 
approximated in the range of to from zero to the dip location (see Fig. [T|) by 7t(to) ps p(m), where we assume 
that J2mP( m ) = 1 as follows from Eq. So the value of the function 7t(to) at the dip is roughly equal to 
p(to*/2) which gives the corresponding maximum of the effective potential V* = — lnp(TO*/2). Thus using the 
Arrhenius law one can expect that the time needed for a random walk to go over this maximum is of order 
r ~ e * = 1/p(to*/2). For p(m) asymptotically behaving as ~ to" - this yields r ~ i.e., a power law in the 
system size N, with a coefficient proportional to ApP . This crude argument gives already a good estimate. It 
can be polished if one implements all details of the zero-range dynamics into the mean-field analysis and works 
out the consequences of the detailed balance condition for the transition rates for a particle to hop into or from 
the monitored node or to stay at it [H, [13] • One obtains an expression for a monitored node i for the waiting 
time T l k ^ m which tells us how long it takes to fall from k to m < k particles, where only the first-passage time 
is taken into account: 

m j m 

r^ m = 22 z2 
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For a regular graph the distribution 7ii(m) is identical for each node so the index i can be skipped. The typical 
melting time, defined to be r = T m „—,o, calculated from this equation grows as before as r ~ N@ , but with a 
slightly modified coefficient which is now proportional to Ap^ +1 . It is only a small correction to the previously 
derived result since the dependence on N is the same. 

In contrast to homogeneous systems, much less is known about the zero-range process on inhomogeneous 
networks, where the symmetry of the partition function {1} resulting from identical node degrees is explicitly 
broken. Some attempts have been made for scale-free networks [19] to explore the phase diagram and the 
dynamics of condensate formation, but not the life time of the condensate. 

We will now take advantage of Eq. (TTTj) to estimate the typical time scale for the dynamics of condensation for 
inhomogeneous systems, assuming that (jll[) applies within the scope of the mean-field approximation also when 
7Tj(m) vary from node to node. We will study the effect of inhomogeneity by introducing to a g-regular graph a 
single node with degree Q > q. This type of irregularity with a single node being different from the others is well 
suited to study the effect of symmetry breaking in zero-range processes, which generates a condensate on a single 
node. The parameter a = h\Q/q plays the role of an external field which breaks this symmetry. Let us first 
determine the static properties of such a system. Because the node weight for the Q node, Pq(to) = Q m p(m), 
differs from the weights for regular nodes, p q (m) = q m p(m), by an exponential factor (Q/q) m , it clear that this 
node has a tendency to attract more particles than the others. We hence expect that 7tq(to) increases fast with 
to, while Tr q (m) decreases. The exact form of these distributions depends on the particular form of the weight 
function p(m) in Eq. |T]) but this does not significantly change the generic behavior. The exponential effect 
coming from the factor Q/q > 1 is dominant. In order to simplify the calculations it is therefore convenient to 
assume that the transition rates do not depend on to: u(m) — 1 or equivalently p(m) = 1. Assuming this, one 
can show 20], that the distribution 7Tq(to) develops a maximum for m close to the upper limit M: 

, . (Q\ m (M + N -m-2\ 
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FIG. 1: Plots of the distribution 7r(m) in the condensed phase, for /3 = 5 and p = 1. The critical density is p c = 1/3. 
From top to bottom: N = M = 50, 100, 200, 400. The position of the dip is marked by a filled circle on each curve. The 
vertical line denotes its asymptotic position as N — > oo: m,/2M = 1/3. 

where we have skipped an overall normalization. As an example, we show in Fig. [5] the distribution 7Tq(to) 
of particles at the Q-node. Above the critical density p c = l/(^ — 1), the distribution has a maximum for a 
number of particles to* which linearly grows with the total system size, to* ~ N, so particles condense at the 
Q-node. In the same figure we also show the distribution of particles on a regular node TT q (m), where it falls off 
exponentially: 

We can now calculate the typical time scales tq = ^ for the condensate to disappear for the first time 
from the Q-node or r q — __,q for the same quantity at the regular g-node, applying Eq. (jlip . The general 
expression for tq is quite complicated [2(j. In Fig. [3] we plot this quantity for N = 20 and various M. In the 
thermodynamic limit for fixed density p it simplifies to 

tq ~ exp {N [p\ogp + pa — (1 + p) log(l + p)]} . (14) 

Furthermore, for p ^> 1 we get tq oc e apN which means that here the characteristic "melting time" grows 
exponentially with the system size, while T q is found to grow only linearly. Actually the time needed for the 
condensate to evaporate from the Q-node can be estimated also using the Arrhenius law, as before, if one applies 
it to the effective potential V(m) — — In ttq (to) . As follows from Eq. (fT2"|) . this function has its maximum at 
to = and the value of this maximum, when one normalizes (|12|) . can be estimated from the Arrhenius law: 
t* « e y (°) = 1/7tq(0) ~ e "ApAf g Q j n con trast to the homogeneous case the time grows now exponentially with 
the system size. 

Let us summarize differences between the condensation observed in zero-range processes on a g-regular network 
and an irregular network. In the first case the partition function is symmetric with respect to permutations of 
the occupation numbers {toi, TO2, . . . , to^} so that any distribution of particles is as probable as any other. This 
symmetry is spontaneously broken at the critical point where a single node containing the condensate becomes 
different from the others. The symmetry is reduced to the group of permutations of the remaining N— 1 nodes. 
On an irregular network, on the other hand, the symmetry is explicitly broken. Because of the nature of the 
symmetry breaking, which produces a condensate on a single node, one can expect that already a model with an 
irregularity on a single node is sufficient to capture the main characteristics of this transition. We have studied 
such a minimal irregularity coming from a single node with degree Q > q. The parameter a = \nQ/q plays the 
role of an external field. The situation is very similar as for standard phase transitions. For a > 0, one observes 
a characteristic exponential suppression n q (m) ~ e~ am of the particle distributions on g-regular nodes which 
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FIG. 2: Plots of tvq {rn) (solid lines) and n q (m) (dashed lines) for the almost g-regular graph with one node Q > q. Here 
q = 4, Q = 16, p = 1 and the critical density p c = 1/3. The curves from top to bottom show: N = M = 50, 100, 200, 400. 
The vertical line marks the asymptotic position of m*/M = 2/3. 
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FIG. 3: The typical time scale tq related to the average life time of the condensate, for iV = 20, q — 4, Q — 16 and 
different M. Circles denote experimental points. 



can be compared to the exponential fall-off of the two-point correlation function in standard field theoretical 
models, while for the absence of the external field one observes long range fluctuations: n q (m) ~ p(m) which 
for the interesting class of weights p(m) is of the power-law type: Tr q (m) ~ mT^ . This change of behavior has 
also immediate implications for the dynamics: a typical time scale for the condensate melting for a > grows 
exponentially with the system size and for a. = sub-exponentially, typically as a power. 

A next step is to consider zero-range processes on complex networks with an arbitrary degree distribution. 
One can in particular address the question of self-averaging that is whether a zero-range process on a single 
typical network chosen at random from the given ensemble of networks with the given degree distribution 
behaves in the same way as the corresponding process averaged over many networks from this ensemble. This 
would be the first step towards the investigations of dynamics on networks coupled to network topology. 
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